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1. Introduction 

Coherence phenomena play an important role in many interesting problems of particle 
physics and cosmology. Well known examples include neutrino oscillations [8, 9], elec- 
troweak baryogenesis [10-16], coherent baryogenesis [6,17], leptogencsis [18-21] and out- 
of-equilibrium particle production [22] . Recently, a novel kinetic theory approach which in- 
cludes both flavour- and particle-antiparticle coherence effects was introduced in refs. [1-7]. 
This formalism, based on the coherent quasiparticle approximation (cQPA) , is in principle 
applicable to all the problems mentioned above. Other recent works related to quantum 
coherence and decoherence in quantum field theory include e.g. [16,23]. 

cQPA scheme is based on the Schwinger-Keldysh formalism of nonequilibrium quan- 
tum field theory [24] with the approximations of weak interactions and slowly varying 
background field {e.g. expansion of the Universe). It provides an extension to standard ki- 
netic theory with quantum Boltzmann equations (see e.g. refs. [14,25,26]) in that in cQPA 
the 2-point correlators are not assumed to be nearly time-translation invariant, which al- 
lows new oscillatory solutions mediating non-local quantum coherence [1-3]. Because of the 
rapid oscillations at quantum scales ~ A;, the collision integrals involving these coherence 
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solutions need to be resummed to all orders of the gradient expansion in Wigner represen- 
tation [4-6] . This resummation procedure ultimately leads to an extended set of Feynman 
rules for coherent systems, as well as coherent quantum Boltzmann equations for enlarged 
set of distribution functions, formulated in [6] for single fermionic and scalar fields. The 
generalization of the formalism to the systems with multiple mixing fermionic and scalars 
fields was presented in [7]. 

In this work, we present a new, more elegant and general, flavour covariant formulation 
of the fundamental cQPA Wightman functions and the perturbative Feynman rules for the 
theory. The cQPA Wightman function is found to be a composite operator, schematically 
of the form iS^^^^ = AaF^^^'^ {t,t)Ab, where the spectral functions A project S^^^ to the 
mixing states labeled by a and b, and the local 2-point vertex function F^i,'^{t, t) contains all 
associated statistical information. The composite structure of the (statistical) propagator 
was also observed in rcfs. [27, 28] for nonmixing scalar and in [21] for fermionic fields 
coupled to thermal bath. In these works, the local correlators F^'^ were enscribing the 
initial conditions of the system, whereas in our approach F^ij'^{t,t) are the dynamical 
variables of interest. This (derived) form can be viewed as a generalization to coherent 
systems of the well known Kadanoff-Baym ansatz [29], which projects the propagator to a 
single mass shell: iSad^ ~ 2Aana, where Ua is particle number in the state a. 

We generalize our previous work also by including the dispersive (hermitian) part of the 
self-energy in the dispersion relations, which extends our formalism to treatment of quasi- 
excitations in plasmas. This is important for example for a consistent approach to neutrino 
oscillations phenomena and thermal resonant leptogenesis [19], as well as for the electroweak 
baryogenesis. We note that the final flavoured kinetic equations for local 2-point functions 
S^''' {t, t), resulting from our approach, arc very reminiscent of the equations of motion for 
the (flavoured) density matrix in the standard density matrix approach (see e.g. [8]). We 
want to emphasize, however, that in our approach all the interaction terms are derived from 
first principles of nonequilibrium QFT. Moreover, these equations are local in time- variable 
t, and they can be written as flavoured quantum Boltzmann equations for the (enlarged) 
set of on-shell distribution functions /(k, t). 

The paper is organized as follows. In section 2, we briefly review the Kadanoff-Baym 
equations and the cQPA scheme. In section 3, we present the flavour covariant derivation of 
the cQPA propagators, and write down the kinetic equations in section 4. In section 5, we 
show the covariant formulation of the cQPA Feynman rules, and in section 6 we use these 
rules to write down the final kinetic equations for the local correlators. In this section, we 
also present a parametric estimate for the error resulting from our approximations. Finally, 
section 7 contains our conclusions. 

2. KB-equations and cQPA 

We consider a spatially homogeneous and isotropic flavour-mixing fermionic system with 
a complex and possibly time-dependent mass matrix m. The complete Kadanoff-Baym 
(KB) equations [29] for the 2-point correlation functions in two-time representation can be 
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written as (cf. [2]) 



-J dt'j:'^''{ti,t')s'^'''{t',t2) = s{ti-t2), (2.1) 

{ij'^dt, - k • 7 - nihih) - ir"ma{ti))S<'>{h,t2) 

- J dif{j:^ih,t')s<'>it',t2) 

+ ^<'>{tut')Stit',t2)) = ±C^°'\h,t2) . (2.2) 

where nih = {m + m))/2 and rua = {m — m^)/{2i), S^'"' are the standard (retarted and 
advanced) pole propagators, iS'"{u,v) = {tp{v)'ip{u)) and iS^{u,v) = {■iIj{u)'iI>{v)) arc the 
Wightman functions, and = {S^ + S"')/2, with similar definitions for the corresponding 
self-energies T,.^ The ± sign in Eq. (2.2) refers to < or > component equations, respectively. 
The collision term is given by 

Cr(ii,i2) = IJ dt'{i:>ih,t')S<{if,t2) - J:<{h,lf)S>{t',t2)) , (2.3) 

and the correlators and self-energies in two-time representation are defined as 

FUti,t2) = y"d=^(xi-X2)e-*("i-''=)F(xi,X2). (2.4) 

Wc will later on refer to Wigner representation as well, where the 2-point functions are 
defined as 

F{k, t) = j d(ti - t2) e*'^°(*i-*^)Fk(ti, t2) , (2.5) 

where k = (A;o,k) is the 4-momentum variable and t = {ti + t2) (2 \s the average time- 
coordinate. 

In cQPA we first expand the KB-equations (2.1-2.2) to the zeroth order in time- 
gradients dtm{t) and in the scattering width V = i{T?' + S^)/2 to obtain flavour-coherent 
propagators with (enlarged) singular quasiparticle phase-space structure. We also assume 
that the hermitian part of the self-energy, Sj^, does not depend on the dynamical correla- 
tors, and thus 9(S|^(ti, ^2) ~ dtm{t) can be neglected as well. The average-time dependence 
of the (dynamical) cQPA-propagators 6*^'^ is found to be factorized to local correlators 
S^'^(t,t) or equivalently to the set of on-shell distribution functions f(k,t). To obtain 
the (first order) equations of motion for S^'^{t,t), we then insert the dynamical cQPA- 
propagators back to Eq. (2.2) and expand to (combined) first order^ in dtm, F and S'*. 
These equations are local in time-variable t, and they can be written as flavoured quantum 
Boltzmann equations for the distribution functions /(k, i). 

'^Note that our definitions of and differs by sign from the more standard convention. 
^Thc dimcnsionlcss expansion parameters are roughly dtm/oj^, Fk/wk and Sk/wk- See section 6.2 below 
for further discussion. 
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3. Structure of propagators 



3.1 Pole propagators and the spectral function 

To the zeroth order in dtm, the non-dynamical pole propagators S""'" are time-translation 
invariant and the KB-equations (2.1) in Wigner representation read^ 

( 1^ -ruh - ij^rua - Tl'{k) ± iV{k))S'''''{k) = 1 , (3.1) 

with a simple solution 

^ ^ ^ — nih — ij^ma — T;^{k) ziz iT{k) ^ ^ ^ 

where we have used S*"'" = T,^ ^ iT and wc keep infinitesimal T. Even though the time- 
gradients dtm and dtTi^{k) are neglected here, m and S^(fc) are understood to depend 
adiabatically on the (average-)time coordinate. Note that the masses mh,a and the self- 
energies S'^ and r carry flavour indices, and thus the computation of the inverse matrices 
in Eq. (3.2) can in general be tedious. The spectral function is then formally given by (we 
now take the limit F — >■ 0) 

A{k) = ^ {S'{k) - S'^ik)) = 7T sgn{ko)6{ p -rrih - il^'ma - ^\k)) 

= 27r '^Aik 6 {ko - ojik) , (3.3) 

i 

where Ai\^ are matrices in Dirac and flavour indices, and Wjk are the energies of the quasi- 
particle states, labeled by index i. We note that in two-time representation the spectral 
function (3.3) satisfies the zeroth order equation 

{ij^'dy -k--f-mh- h''ma)AUy) - J dy' S(^(y - y') A(y') = . (3.4) 

where y = ti — 1,2 is the relative time-coordinate. Moreover, from the singularity structure 
of Eq. (3.3) it clearly follows that Ak{y) is analytic in y, and we can therefore expand the 
convolution integral in Eq. (3.4) as 

/,. oo 
dy' ^^y - y')AUy') = / dy' ^^{y')AUy -y') = Yl ^kn ^''^;AUy) , (3.5) 
n=0 

where 

St^^y dyy^^^iy). (3.6) 

Using this expansion in Eq. (3.4), and solving iteratively for dyAk{y) neglecting the con- 
tributions of order (E^)^ and higher,^ we find 

dyA^iy) = ~iH^,esAUy) => A^{y) = e-'^^^'^'^v A^,o , (3.7) 



^If the average-time derivatives dtS^'°'{k, t) would have been kept in Eq. (3.1), they would have eventually 
been found to vanish in the limit of vanishing dtm and dt'E{k). 

*This is consistent with the general procedure of QPA, where the fco-poles of the spectral function (3.3) 
are solved iteratively by using the lowest order solution for ko inside T,'^{ko, k). 
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where we denote Ako = .4k (0) and the effective quasiparticle Hamiltonian is given by 

^^k,eff = H^o + J dy 7°Si^(y) e^^"-"^ , (3.8) 

where -ffko = k • a + 7*^m/j + i^^^^ma is the free Hamiltonian of the system. Note that 
we have not used the spectral sum-rule: .4k(0) = J ^A{k) = ^'J^Sij, here, because in 
general it does not hold exactly in QPA with pole contributions only [30] . It is now trivial 
to show that the quasiparticle spectral function satisfying Eq. (3.4) has an important time- 
additivity property: 

Ai,{yi + y2)A^o = AUvM^o^MA^o , (3.9) 

such that 

uUy) ^ My)A^o = e"'^"-"^ (3.10) 

can be interpreted as a time-evolution operator. Note that Uk{y) is not necessarily unitary 
in quasiparticle approximation due to possible breaking of the spectral sum-rule. How- 
ever, if the sum-rule holds, it follows from the general hermiticity property: Al^{y) = 
7°-4k(— y)7'^, that i?k,eff must be hermitian and the "free" (quasiparticle) evolution is uni- 
tary. We want to emphasize that the additivity property (3.9) holds irrespective of the 
unitarity of U]^{y). This property is crucial for showing below that the dependence on the 
relative coordinate y = ti—t2 and and the average coordinate t = (ti + t2)/2 of the Wight- 
man functions 5^'^ factorizes neatly in quasiparticle approximation, leading to (enlarged) 
singular phase-space structure and local time-evolution equations for the corresponding 
on-shell distribution functions. 

3.2 Wightman functions 5^'^ 

To the zeroth order in dtm and T ~ 5]^'^, the KB-equation (2.2) for the Wightman 
functions: 

(i7% -k-7-m/,-i75m„)5<'>(ti,i2)- J dy' J^^it^ - t')S<'> (t' ,t2) = , (3.11) 

is identical to Eq. (3.4) for the spectral function with the replacement y ^ ti, and t2 

acting as a parameter. Therefore, it has a factorized general solution (matrix product): 
5^'^(ti,t2) = A]i{ti)C^'^ {t2), where t2-dependence can be fixed by using the hermiticity 
properties: i5<'>^(ti, is) = 7°«5k '^(*2, ii)7° and Aliy) = 7°-4k(-y)7°, to get 

S^'^iti,t2) = Ai,{ti)A^^S<'> (0,0) A^^A^,{-t2) . (3.12) 

Eq. (3.12) combines the zeroth-order time evolution in the relative coordinate y = t\ — t2 
and in the average coordinate t = (ti + i2)/2. This general structure of the statistical 
propagator was also observed in [27,28] for nonmixing scalar and in [21] for fermionic field 
coupled to thermal bath. 

To factorize y- and t-dependencies in Eq. (3.12), we use the additivity property (3.9) 
and Eq. (3.12) for local correlators S^'^{t,t) to get 

S^'^iti,t2) = Ai,{y/2)A^^S<'>{t,t)A^^Ai,iy/2) . (3.13) 
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In Wigner representation Eq. (3.13) becomes a convolution integral: 

dk' 



/dk 
-^A{k'o, k)^-(}5<'>(i, t)A^^Ai2ko - k'o, k) 

= 2TrJ2 ^ik^M '^(*' t)A^o^jkS(ko - ^{ojik + cojk)) , (3.14) 



where we have used Eq. (3.3) for the spectral function. This structure of cQPA Wightman 
functions is one of the main results of this paper. The local correlators S^''^{t,t) con- 
tain the full statistical information of the system, including phase-space particle number 
and (flavour) coherence distribution functions. Wc note that if A^QA{k) commutes with 
A^qS^'^ {t,t) for all ko, the time-additivity property (3.9) yields S^'^{k,t) = 
A{k)A^QS^'^ {t,t). Moreover, in the absence of the dispersive self-energy S'* and the co- 
herence contributions (see Eqs. (3.16-3.17) below), we recover the standard form of the 
Kadanoff-Baym ansatz: iS^''' {k,t) oc 2/^''' {k,t)Aii{k), where f^'^{k,t) arc the phase- 
space distribution functions of flavour state i. Based on these considerations, Eq. (3.14) 
can be viewed as a generalization of the Kadanoff-Baym ansatz to systems including flavour- 
and particle-antiparticle coherence. 

We also note that the composite structure of the Wightman function (3.14) involves 
cross-couplings between different flavour-components of 5"^'^(A;, t) and S^'^{t,t), if the 
spectral function is not flavour diagonal.^ For example, in the case of two-flavour mixing 
we have: 

^ ^(^(fc[„k)^^i)^^5<;>(i,i)Ko'^(2fco - 

+ {Aik'„k)A^^)^^S<;{[{t,t){A^^A{2ko - k'o,k))^^ 
+ {A{k'„k)A^^)^^S<;>it,t){A^^A(,2ko - k'„k)). 

and similarly for other components S^'^{k,t). 



(3.15) 



3.3 The limit of vanishing dispersive self-energy 

In the limit of vanishing dispersive self-energy E'*, the spectral function (3.3) is diagonal 
in the mass-diagonal basis (see [7] for details): 

Aij{k) = TT sga{ko)0 + mi)6{k^ - mj)5ij , (3.16) 

where rrii are (positive and real) elements of the diagonalized mass matrix ma = UmV^ , 

U and V arc unitary N x A'"- matrices, and the mass-basis correlators and sclf-cncrgics arc 
defined as Sd{k,t) = YS{k,t)X^ and Sd(/c,t) = XS(A;,t)yt with Y = Pl®U + Pr®V 
and X = 7*^1^7*'. Using the spectral function (3.16) in Eq. (3.14) and parametrizing the 
most general (spatially homogeneous) local correlator <S'^'^(i, t) in terms of the helicity and 



^Note that the cross-couplings always vanish in fco-integration, reducing Eq. (3.14) to trivial identity for 
local correlators. 
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energy projection matrices: -Pfe(k) = |(l + h{k ■ 7°77^)/|k|) where h = ±1, and 1^.^ + rui, 
we get 



- + mj)f^;f:±S{ko T Aa;,,)l , (3.17) 



where kf^ = {±0Ji,k), a;i(k) = + mf , coij = {oji + ojj)/2 and Aojij = {coi - coj)/2. 

The Wightman functions S^'^{k,t) have a singular phase-space structure with the shells 

ko = ±uJij and k^ = zizAcoij multiplied by the corresponding projection matrices and the 
on-shell distribution functions f^j^^ (k,t) and f^^'^ (kjt), which describe the amount of 
flavour coherence between the flavour eigenstates with energies ibwj and ibcjj, respectively. 
For more detailed discussion and an alternative derivation of the Wightman functions 
(3.17), see [7]. 



4. Kinetic equations I 

The first order KB-equations (2.2), written for the hermitian local correlators S^'^(t,t) = 
zS'^'^(i,t)7° are given by^ 

dtS<'>{t,t) = -i[H^o{t),S<'>it,t)]-i J dt' {tUt,t')S<'>{t',t)- S<'>{t,t')tUt',t)) 

±7°(C£°^i,*)+Cr(i,i)^)7°, (4.1) 

where S'* = 7'^S'*. Here we have dropped the interaction terms ~ S'^'^S'^, which arc 
related to finite scattering widths of the propagators and arc therefore beyond quasiparticle 
approximation. The kinetic equations (4.1) arc non-local in t, because the interaction terms 
(convolutions) involve non-local correlators Sl^{t',t). Moreover, the loop expansion of the 
self-energies S involves in general (non-local) correlators Sj^{wo,WQ), where vertex time- 
arguments Wo and w'q are integrated over. However, local kinetic equations can be obtained 
up to first order by using the zeroth order equations (following from Eqs. (3.12) and (3.9)): 

5<'>(w;o, w'o) = AUwo - t)A^^S<'>{t, t)A^^AUt - w'o) , (4.2) 

to expand all non-local correlators in terms of the local correlators S^'^{t, t) to get a closure 
for the kinetic equations (4.1). This expansion of the interaction terms can neatly be done 
in momentum space by enlarged set of Feynman rules involving flavour-coherent effective 
propagators. 

The hermitian conjugate terms arise from the Umit dtS^'>{t,t) = Wmt^^t, {dt^S^'> {ti,t2) + 
dt2S^'^ {t\,t2)) = \vcat2^ti dti{S^'^ (t-i,t2) + t2)) , or alternatively Eq. (4.1) can be obtained 

as fco-integral of the antihermitian part of KB-equation (2.2) in Wigner representation. 
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5. Momentum space Feynman rules 



To derive the cQPA Feynman rules, we first rewrite the interaction convolutions in Eq. (4.1) 
to the form 

iit)= J dt'i:ut,t')s<'>{t',t) 

= J dwodw'odw'^ 5{t - w;o)Sk(w;o, w'o)S^'>{w'o, w^w'^ - t) . (5.1) 

Then we expand all Wightman functions in terms of S^'^{t, t) using Eq. (4.2), and perform 
a shift of integration variable: wq — > wq = wq — t, for the time-arguments wqjw'q, Wq and all 
vertex time- arguments inside Ti{wo,Wq). As a result all Wightman functions in Eq. (5.1) 
are replaced by (effective) propagators 

S^'^iwo,w'o) = S^'''{wo + t,w'o + t)=Ak{wo)A^^S<'>it,t)A^^Aki-w'o), (5.2) 

while the translation invariant pole-propagators appearing in the generic loop diagrams in 
Closed Time Path (CTP)-formalism [24] remain invariant, since S^'"' {wq — Wq) = S^'°'{wo — 
w'q). The primary CTP-propagators, given by^ 

g++ = st= s'' -s<, s+- = -S< , 

S— = 5* = -S" - S< , S-+ = S> , (5.3) 

then transform accordingly. Similar transformation properties apply to (coherent or trans- 
lation invariant) propagators of the other fields coupled to the field ip, appearing inside 
the self-energies. By writing the effective propagators (5.2) as double Fourier-transforms 
w.r.t wq and w'q, and performing the time integrals over all vertices wo,w'q, w'q and the ones 
inside T,k{wo, w'q) = Sk('i(;o + t,w'Q + t), we then obtain 

I{t)= J du)od4d4'5(u;o)Sk(u'o,4)'5^'^(4,4')^K) 

= -| ^iSeff(fe,tXffi(fc,t), (5.4) 

where the effective momentum space self-energies are formally given by^ 

^es{k,t)= J dw'Qe-'''<^<t^{0,w'Q) 

= J dw'Qe"''>^'-<^^k{t,w'Q), (5.5) 

while the effective in-propagators are defined as 

iS^s%{k,t) = A{k)A^^iS<^>{t,t) . (5.6) 

''identical relations hold for the self-energy components E"**. 

®Note that this relation is not the standard Wigner transform of the self energy. 
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The effective Wightman functions inside T,es{k,t) are written as 

iS^4>{k, k', t) = A{k)F^^{t, t)A{k') , (5.7) 

with 

F<^>it,t) = A^^zS<'>{t,t)A^^ (2vr)3<53(k - k') . (5.8) 

Moreover, S^^f' {k , k' , t) are accompanied by an additional 4-momcntum integral w.r.t k' . 
On the other hand, the pole propagators inside Seff(fc, i) have the standard expressions, 
given by Eq. (3.2). The effective self-energies Ties{k,t) are then simply the usual (CTP) 
self-energies in momentum space, where all propagators are replaced by the corresponding 
effective propagators. We note that Eq. (5.4) is of the same form as the (standard) inter- 
action term arising from zeroth order truncation of the (O-)gradient expansion in Wigner 
representation (see ref. [7] and Eqs. (6.12-6.13) below). However, Eq. (5.4) with effective 
propagators and self-energies involves a resummation of rapid coherence oscillations of fre- 
quencies ~ k, contributing to all orders in the 0-expansion. Only in the case of vanishing 
coherence contributions, Yles{k,t) and (flavour diagonal) S^g^^{k,t) reduce to standard 
Wigner representation expressions 'E{k,t) and S^'^{k,t), and I{t) reduces to the zeroth 
order truncation of the 0-expansion. 

Furthermore, the hermitian conjugate interaction terms in Eq. (4.1) are given by^ 

= J ^r^5e1^,out(^' th^J^Uk, t) , (5.9) 
where the effective out-propagators are given by 

iS<^>^,ik,t) = -7°<ff2(fc,t)7° = iS<'>{t,t)A^^A{k) . (5.10) 

To summarize, wc list the cQPA Feynman rules to compute the interaction terms I(t) and 
I^{t) of Eq. (4.1) in momentum space in terms of (effective) flavour-coherent propagators 
involving the (dynamical) local correlators S^'^{t,t): (see Fig. 1) 

• Draw the self-energy diagram(s) for zSefr(A;, f) and use the standard CTP-Feynman 

rules for the vertices in momentum space, and apply CTP-indices for the propaga- 
tors. The momentum conservation delta function is not applied to the out- vertex of 

iSeff(fc,t).10 

• Use the relations (5.3) to write the CTP-propagators in terms of iS'^'"' and iS^'^ . 

• For the pole propagators, use the standard expressions iS'^'°-{k) of Eq. (3.2) (with 
r ^ 0) with the usual (single) /c-integration. 

• For the internal Wightman propagators, use iS^4>{k,k',t) of Eq. (5.7) with double 
4- momentum integration over k and k' . 

• For the external in- and out- Wightman propagators, use iS^^^^ik, t) and iS^Q^^^^{k, t) 
of Eqs. (5.6) and (5.10), respectively. 
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Internal propagators: 

iS^s'^ik, k\ t) = A{k) F<^r{t, t) A{k') 

In- propagators: 

^S<s%{k,t) = A{k)A^,]^S<^>{t,t) 

Out-propagat ors : 

^S:s>{k,t)=^S<^>{t,t)A^^A{k) 



Figure 1: The Feynman rules for the flavour-coherent (effective) Wightman propagators. The 
propagators are composed of the spectral function A{k) (fermionic lines) and the effective 2-point 
vertices (circled crosses) encoding the quantum statistical information. 



0< ® (B^^i 
'^Ses;m{k,t) ^'S'gfj'o^t(A;, t) 

^^es{k,t) 70zElg.(A;,t)70 

Figure 2: A diagrammatic representation of the interaction terms I{t) (left) and 7*'/^(f)7° (right) 
in the kinetic equation (4.1). 




In figure 1 we show diagrammatic representations of the flavour-coherent cQPA propa- 
gators, which are composed of the spectral function A{k) (3.3) and the effective 2-point 
vertices F^^{t,t) (5.8) (or A^^iS^'^ {t,t) and iS^''^ {t,t)A^Q for in/out-propagators, re- 
spectively), encoding the quantum statistical information. In figure 2 we show generic 
structures of the interaction terms I{t) and I^{t) appearing in the kinetic equation (4.1). 

These Feynman rules can be used to compute any type of Green's function in OTP 
formalism, not only the self-energies appearing in (KB-) kinetic equations. For example, a 



®Note that the effective self energies Ecft (fc,t) do not follow the hermiticity properties of the Wigner 
transformed self-energies E(fc,f). 

^"in comparison, by the standard (momentum space) Feynman rules, any vertex can be singled out to be 
lacking the momentum conservation delta function. 
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perturbative contribution to full Wightman function becomes: 
iS<^{k,t)= j d(ti - t2)e^'=°(*i-*^) j dw;odw;(,i-S<(ti,w;o)zS^(u;o,w;())i5<K,i2) + ... 
=^ / ^(04^5eff((4k),P,t)iS> (p,i)i5< (p,(2A;0^ 

(5.11) 

where we have again used shifts in the integration variables such as wq ^ wq = wq — t. 
We notice the similar convolution structure as in the case of zeroth order correlator in 
Eq. (3.14). For the corresponding local correlator we then get by fco-integration over 
Eq. (5.11): 

iS^,^.At.t) = j ^.5<,„„t(fc,*)^S>ff(fe,t)^5<,in(^,^) + -- - • (5.12) 

The perturbative expansion in Eqs. (5.11-5.12) is performed around the time instant t, 
i.e. the implicit local correlators S^''' {t, t) appearing on the RHS carry the full (quantum) 
statistical information. Therefore Eq. (5.12) is actually an equation for S^(t,t). On the 
other hand, by performing the perturbative expansion around the "initial" local correlators 
£'^'^(0,0), we find a relation 

iS<,,{k, t) = A^,{t)A^^iS<,,{k, 0)A^^AU-t) , (5.13) 

where iSp^j.^{k,0) can be read off Eq (5.11) by setting t = 0. By identifying C/k(t) = 
A]i{y)A^^Q as a quasiparticle time-evolution operator according to Eq. (3.10), we see that 
Eq. (5.13) corresponds to perturbative time evolution in quasiparticle interaction picture. 
It is well known that out-of-equilibrium perturbative time evolution suffers from secularity, 
and the error becomes large for Ft > 1 (see e.g. [31]). In contradistinction, in cQPA we 
do not use the perturbative expansion around initial time correlators. Instead we expand 
the interaction terms I{t) in the kinetic equation (4.1) in terms of the local correlators 
S^'^{t,t), and solve the resulting local time-evolution equations for S^'^{t,t). Since the 
expansion is made around the time-evolved correlators, this procedure does not suffer from 
secular error growth. We will show this in section 6.2 below by computing iteratively the 
magnitude of the neglected correction terms. 

5.1 Example: One-loop self-energy in Yukawa theory 

As an example of the Feynman rules, we write down the effective one-loop self-energies 
iE^^(/c,t) in a Yukawa theory described by the (interaction) Lagrangian 

Ant = -yij i^i (p Phil^j + h.c. , (5-14) 

where yij is a complex Yukawa coupling matrix in flavour indices. The one-loop self-energy 
diagram(s) can be obtained for example from the two-particle irreducible (2PI) effective 
action [32] two-loop vacuum diagram: 

Espi = -yijy\m I IV [PLSjiiu, v)PRSmi{v, u)] A(u, v) , (5.15) 

Jc 
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where sum over repeated flavour indices i, j, I, m is understood, and the integration is over 
the Keldysh path (see e.g. [1]). A{u,v) denotes the scalar field (CTP) propagator. The 
fermion self-energy now follows by a functional differentiation: 



dS°f{v,u) 



=i [yjiymjPRSfiiu, v)PlA'"'{v, u) + v)PrA''\u, v)] . (5.16) 

By performing a Wigner transformation, we now find for the self-energies iT,^'^{k,t) 

iT<>{Kt)= j (2^)4 (2^^4 i^^f6\k -q-p) [yliym.PRiS^r;,^ {q, miA>'<{-p, t) 

+yuyljPLiSf^''iq,t)PRiA<'>{p,t)] . 

(5.17) 

On the other hand, the effective self-energies iS^g^(fc, t) appearing in the kinetic equations 
(4.1) (and below in Eqs. (6.1) and (6.7)) are obtained by using the cQPA Feynman rules de- 
veloped in this section. By assuming that the scalar propagators A^'-^ are non-coherent,^^ 
we get (see Fig. 3) 

X [yliymjPRiS^iiM , q, t)PLiA>'<{-p, t) 
+yuyl,PLiS^dM(^',q,t)PRiA<'>{p,t)] . (5.18) 

We notice that the 4-momentum is not formally conserved in the out-vertex of this self- 
energy diagram (unless q' = q). However, the ^'-integration can always be performed 
trivially because of the singularity structure of S^^^ (q' ,q,t). This results to appropriate 
energy projections of the spinor structure, such that the kinematical interpretation of 
the (coherent) collision process remains natural. By using the explicit expressions (5.7- 
5.8) one can write Eq. (5.18) in terms of the local correlators 5q'^(t, t), and perform 
some of the momentum integrations. Further simplifications of Eq. (5.18) would require a 
specification of the fermionic spectral function (3.3) and the scalar correlators A"^'^ (p,t), 
and are therefore model dependent. 



6. Kinetic equations II 

In this section, we use the cQPA Feynman rules to rewrite and simplify the interaction 
terms in the first order kinetic equations (4.1) in order to obtain (local) time-evolution 
equations for the local correlators S^{t,t). In the case of vanishing dispersive self-energy 
S'*, we also write down the corresponding flavoured quantum Boltzmann equations for the 

^^The single flavour scalar field could also mediate (flavour diagonal) particle-antiparticle coherence (see 
refs. [3,6]). The non-coherent scalar propagators A*~'^(p, t) correspond to the usual quasiparticle propaga- 
tors (see e.g. [30,33]), possibly involving out-of-equilibrium particle number distribution functions. 
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iA>'<(5-fc,<) iA<'>(fc- 

/ \ / \ 

k ' q' q ^ k k ' q' Q ^ k 

< I < (g) < I < < I « ® < I < 

Figure 3: The Feynman diagrams contributing to the (effective) self-energies at one- 

loop level {i,j,l,m are flavour indices). 



on-shell distribution functions, given by Eq. (3.17). Finally, we present an estimation for 
the error in time evolution resulting from our approximations. 

To begin, we use the results (5.4) and (5.9) with Eqs. (5.6) and (5.10) for the interaction 
terms in the kinetic equation (4.1) to get for 5*^ 

dts<{t,t) = -i[H^Mt),s<{t,t)r - 1 {{t>{t),s<it,t)r - {t<{t),s>{t,t)r) , 

(6.1) 

where we have defined 

HkMi) = Hko{t) + f:Ut), (6.2) 
and (for all components of self-energy) 

Sk(i) = ^cs{k, t)Aik)A^^ = ^esi^^n., K i)^ik^ko ' (6-3) 

i 

where we denote T,^ = j^T^^ and S^'^ = ^^iH^'^ and we have used Eq. (3.3) for the 
spectral function. The generalized (anti) commutators (with respect to both Dirac and 
flavour indices) are defined as 

[A, B] = AB- B^A^ , {A, B}"" = AB + B^A^ . (6.4) 

We don't need to consider the kinetic equation for {t, t) , because the Wightman functions 
are in general related by 

S> = -S< - 2iA (6.5) 

and the spectral function A, given by Eq. (3.3), is non- dynamical. The effective self-energies 
Sgff (fc, t) are written in terms of the local correlators S^{t, t) (and similarly for other fields 
coupling to field '0) using the Feynman rules of the previous section. The kinetic equation 
(6.1) is therefore local in time t for arbitrary interactions, and hence an ordinary (coupled) 
first order differential equation in contrast to integro-differential equations with full memory 
integrals (see e.g. [34]). However, the kinetic equation (6.1) is not in general diagonal in 
momentum |k| due to (possible) self-interactions. Furthermore, if T,^ does not depend on 
the dynamical correlators, we find to leading order in S'*: S|^(t) = J dy E^(y, t)e*^''0^*)^ 
such that -ffk,eff(0 = -f^k,eff(^)) where the adiabatic t-dependence is now explicitly shown. 
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We have assumed this non-dynamical S generahy in the zeroth order equations in order 
to solve the phase-space structure of the cQPA propagators. However, in the (first order) 
kinetic equation (6.1) we can relax this assumption and consider dispersive self-energies 
S'* which involve dynamical correlators and (flavour) coherence in the same footing as 
absorptive self-energies S^'^. 

In the case of an interaction with a thermal bath, the leading order self-energies do 
not depend on the dynamical correlators (thus T,es{k,t) — >■ T,{k,t)) and they satisfy KMS 
relations [35] (/3 = l/T): 

J]>{k,t) = e^^°E<{k,t). (6.6) 
Using the relations (6.5-6.6) in the kinetic equation (6.1), wc then find a linear equation 

at5<(t,t)=-z[iJk,efr(o,^k -({fk(i),^^(i,t)r- / ^{f{k,t),s<^ik,t)r), 

(6.7) 

where T-kit) is given by the analogous relation to Eq. (6.3) with T = 7°r = (S> -|- E<)/2, 
and 

iS<{k,t) =2f^^{ko)A{k,t) 

=277 sgn{ko)feq{ko)d{ -ruhit) - ij^ma{t) - ^^{k, t)) , (6.8) 

with feq{ko) = (e^^° -|- 1)~^, is the standard quasiparticle Wightman function in (local) 
thermal equilibrium [33]. 

The structure of the kinetic equation (6.1) or (6.7) with generalized (anti)commutators 
is very similar to the equations of motion for the flavoured density matrix in the standard 

density matrix approach to flavour mixing phenomena (see e.g. ref. [8]). Indeed, the local 
correlator S^-j{t,t) is closely related to the density matrix in the basis of eigenstates 
ti*(k, /i)a|i^|0) and Vi(k, h)b^-^^\0). However, to obtain the usual equation of motion for the 
integrated density matrix in the flavour indices only, one would need to make a factorization 
ansatz in Eq. (6.1): 

S^,ij{^,t)^9{^)Pij{t), (6.9) 

where the dependence on the momentum k and Dirac-indices is factorized to the function 
^(k). In general, this kind of ansatz is not justified unless some stringent extra constraints 
are imposed. Wc want to emphasize that our kinetic equations (6.1) and (6.7) with the self- 
energy functions (6.3) are derived from first principles using the methods of nonequilibrium 

quantum field theory. 

In order to solve the kinetic equation (6.1) for S^{t,t), one needs to parametrize the 
most general spatially homogeneous local correlator S^{t, t) in terms of eight independent 
matrices of (sub-)Dirac algebra, consisting of the products of 7°, k • 7 and 7^. Then, 
using the hermiticity properties (under the Dirac and flavour indices) wc obtain coupled 
equations of motion for 8A^^ real functions of |k| and t, where N is the number of flavours. 
We note that the equation (6.1) can be solved in any flavour basis, where the quasiparticle 
spectral function (3.3) is determined. 



- 14 - 



6.1 Flavoured quantum Boltzmann equations 

In the limit of vanishing dispersive self-energy S'*, we can use Eq. (3.17) for the Wight- 
man functions in the mass basis, where the on-shell distribution functions /^^^^(k, t) and 
fijh±0^i^) the dynamical variables. Using the projection operators Ph and + rrii, 
we can then solve the distribution functions in terms of local correlators S^''^{t, t) and use 
the kinetic equation (6.1) to obtain 

dtf^^^ = T i2Aa;.,/™<^> + + C±[/]., , (6.10) 

^*/5± = T ^2cD,,-/S±^ + iX'h±[f]ij + Cf,Af]ij , (6.11) 

where X^_^[f] involve terms proportional to mass and mixing gradients: dtm.i/oj'j and 
'E'^ = i{ydtV^ ± UdtU^)/2, and C^'^f/] are the collision integrals. Prom these flavoured 
quantum Boltzmann equations (6.10-6.11) we can directly see that to zeroth order the on- 
shell functions flj^^ and f^^'^ are oscillating at frequencies Autij and Oij, respectively. 
We present the full expressions of X^_f_ and C^_|_ in [7] for a general N x A^-mixing scenario 
and a general Dirac decomposition of spatially homogeneous self-energies. 

6.2 Estimation of error in time evolution 

In order to obtain the local evolution equation (6.1) for S^''' {t, t) we have extensively used 
the zeroth order equations (3.9) and (3.12) inside the interaction terms I(t) of the full 
kinetic equation (4.1). In this section we estimate the magnitude of the error resulting 
from this approximation. 

First, we notice that the convolutions I{t) can be written in Wigner representation as 

/fuii(i)= / dt'^i,{t,t')S<'>it\t) 

= J ^e-^0{nk,t)}{S<'>ik,t)} 

= I ^i;^(|)>((5^"oS(^>o)5<'>(fc,t)) , (6.12) 

n=0 

where the 0-operator is defined as 

0{/}{<7} = \{dtfdk,g - dkjdtg) . (6.13) 

Next, by using Eqs. (3.14) and (4.1) or alternatively KB-equation (2.2) in Wigner rep- 
resentation, we can parametrically estimate the time-derivatives appearing in Eq. (6.12) 
as 

dtS{k, t) ~ koS{k, t) + e koS{k, t) , 

dtT.{k, t) ~ feoS(fe, t) + e koT,{k, t) , (6.14) 

with 

(dtm S(M)\ 
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The zeroth (e-)order terms in the RHS of Eqs. (6.14) arise from the commutators of the 
(coherent) correlators 5^'^ with the Hamiltonian ifko in Eq. (4.1), and accordingly for 
the self-energies involving coherent correlators. We may further estimate dkQT,{k,t) ~ 
S(A;, t)/ko to get 



where we have noticed that all zeroth order terms have been accounted for in the effective 
propagators and self-energies. We see that the corrections to the leading order terms in- 
cluded in kinetic equations (6.1) are suppressed by powers of ^jm/feg and S/Zcq. Therefore, 
we conclude that in the case of weak interactions and small gradients of the background 
field (e <^ 1) the error of expanding the interaction terms I{t) in terms of local correlators 
S^'^{t,t) by using the cQPA Feynman rules is small for any instant of time t. Conse- 
quently, the error of time evolution is not growing as Tt as in the case of secular evolution, 
but rather, our error is that the scattering rate F (and the thermalization time ^ 1/r) is 
misevaluated by a small correction: Ff^n « (1 -|- e)F. 

7. Conclusions and outlook 

In this work we have presented a flavour-covariant formulation of the coherent propaga- 
tors and Feynman rules for flavour-mixing fermionic fields in spatially homogeneous and 
isotropic systems. Our formalism is based on the coherent quasiparticle approximation 
(cQPA) [1-7], where nonlocal coherence information is encoded in new spectral solutions 
at off-shell momenta. In addition to more simplified covariant formulation, we have gener- 
alized our previous work [7] to the case of nonzero dispersive self-energy, which leads to a 
broader range of applications. 

The new formulation is based on the two-time representation of the correlators and 
it immediately reveals the composite nature of the cQPA Wightman function as a prod- 
uct of two spectral functions and a local correlator, interpreted as an effective two-point 
interaction vertex, which contains all quantum statistical and coherence information (see 
Eqs. (3.13-3.14) and (5.7)). This derived form can be viewed as a generalization of the 
well known Kadanoff-Baym ansatz to coherent systems. In section 5, we have derived the 
cQPA Feynman rules involving the flavour- coherent (composite) propagators, presented in 
Fig. 1. These rules are designed for computing the interaction terms of the Kadanoff-Baym 
equations (4.1) in momentum space in terms of the local correlators S^''^{t,t), and they 
can in principle be applied to arbitrarily complex self-energies involving flavour-coherent 
propagators inside the loop-integrals. 

We have used the cQPA Feynman rules to derive flavoured kinetic equations (6.1) for 
the local correlators S^'^{t,t) in the presence of dispersive and absorptive interactions. 




(6.16) 



and furthermore by using Eq. (6.12) 




n=l 



(6.17) 
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These equations are local in time- variable t, and they can as well be rewritten as flavoured 
quantum Boltzmann equations for the (enlarged) set of on-shell distribution functions 
/(k, i). The structure of Eqs. (6.1) is very intuitive, consisting of the generalized commu- 
tators and anticommutators in both flavour and Dirac indices. They are reminiscent of the 
equations of motion for the (flavoured) density matrix used e.g. in refs. [8]. We want to 
emphasize that in our approach all the interaction terms, given by Eq. (6.3), are derived 
from first principles of nonequilibrium quantum field theory by consistent approximations 
using the cQPA scheme. We have further shown that the corrections to our results are 
(parametrically) suppressed by the powers of dtm/kQ and T,{k,t)/ko. 

In the present work we have extended the cQPA scheme to include dispersive self- 
energy corrections, which allows the treatment of quasi-excitations in plasmas. The natural 
applications of the formalism include for example neutrino flavour oscillations in the early 
universe [36] and clcctrowcak baryogenesis, where the formalism needs to be generalized 
to the case of a stationary planar symmetry, as discussed in ref. [1]. Another interesting 
application to our formalism is a consistent approach to resonant leptogenesis, where the 
flavour coherence effects of the singlet neutrinos are expected to play an important role in 
the dynamics. In the regime of nearly degenerate neutrino masses, the formalism should 
be extended to include appropriate finite width corrections for the propagators. 
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